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We study a ID system with a power-law quasiparticle dispersion oc |fc|“ sgnfc in the presence of 
a short-range-correlated random potential and demonstrate that for a < 1/2 it exhibits a disorder- 
driven quantum phase transition with the critical properties similar to those of the localisation 
transition near the edge of the band of a semiconductor in high dimensions, studied in Refs. [I] and 
[2]. Despite the absence of localisation in the considered ID system, the disorder-driven transition 
manifests itself, for example, in a critical form of the disorder-averaged density of states. We 
confirm the existence of the transition by numerical simulations and find the critical exponents and 
the critical disorder strength as a function of a. The proposed system thus presents a convenient 
platform for numerical studies of the recently predicted unconventional high-dimensional localisation 
effects and has the potential for experimental realisations in chains of ultracold atoms in optical 
traps. 
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It is generally believed that increasing disorder 
strength in a conducting material in dimensions d > 2 
leads to the Anderson localisation transition^ with uni¬ 
versal properties that depend only on the space dimen¬ 
sionality d. 

However, as we have demonstrated recently[TJ 2], in 
high dimensions d > d c (with d c = 4 for conventional 
weakly doped semiconductors and d c = 2 for Dirac 
semimetals) the phenomenology is significantly richer. 
Namely, a material with a power-law quasiparticle spec¬ 
trum oc k a in the presence of a short-range random 
potential in high dimensions d > 2a exhibits an un¬ 
conventional disorder-driven quantum phase transition in 
the bottom of the band, that lies in a universality class 
distinct from the Anderson transition jT] [2j. Almost 30 
years ago|!H6j the existence of such a transition was sug¬ 
gested for the specific case a = 1, d = 3 for 3D Dirac 
materials that have been later extensively studied in the 
literature [U |2J ITHTT] establishing a consensus for the ex¬ 
istence of this novel transition in 3D Dirac semimetals. 
Recently [2 [2] we have shown the existence of such an 
unconventional disorder-driven transition and studied its 
properties for arbitrary a and d, such that d > 2a, 
demonstrating that it is a generic property of high di¬ 
mensions and is not specific just to Dirac semimetals. 

For materials in d > 2 dimensions in the symmetry 
classes that allow for localisation, this transition coin¬ 
cides with the localisation transition for the states in 
the bottom of the bandjSJ- However, it exists even if 
all states are always localised (e.g., in d < 2 dimen¬ 
sions) or if localisation is disallowed by symmetry [e.g., 
in 3D Weyl semimetals (WSMs) with sufficiently smooth 
disorder!!!), [T3]] and manifests itself, for example, in the 
critical behaviour of the density of states and conductiv¬ 
ity. 

So far such unconventional disorder-driven transition 


has yet to be observed experimentally. Perhaps the main 
obstacle in 3D Dirac materials is the long-range Coulomb 
nature of quenched disorder distinct from short-range 
random potential required to observe the critical be¬ 
haviour of the conductivity pLj (although the transition 
in the density of states is still observable in the pres¬ 
ence of Coulomb impurities). Another possible plat¬ 
form for studying high-dimensional localisation phenom¬ 
ena is periodically-kicked quantum-rotor systems, that 
can be mapped [21 E] onto high-dimensional semicon¬ 
ductors with quadratic spectra. Such systems have been 
used to simulate 1D{I6], 2D[17], and 3D[miig Ander¬ 
son localisation, but the case of higher dimensions still 
remains to be realised. 

Numerical simulations in high dimensions may be ex¬ 
tremely demanding in terms of computing power. For in¬ 
stance, the quadratic spectrum of long-wave excitations, 
generic for lattice models with short-range hopping and 
inversion symmetry, corresponds to d c = 4 and thus re¬ 
quires simulations in d > 5 dimensions, with the number 
of sites growing rapidly oc L d as a function of the linear 
size L of the system. 

In this paper we suggest and study a new playground 
for unconventional disorder-driven transitions, which is 
rather convenient for numerical simulations and is also 
currently accessible for experiments: ID systems with 
long-range hopping. 

Because the concept of high dimensions d is defined (TJ 
[2] relative to the quasiparticle spectrum, the physics of 
high-dimensional disorder-driven transitions can be ob¬ 
served in any dimension d by appropriately designing 
the inter-site hopping; e.g., realising the quasiparticle 
spectra oc |fc|“ with a < d/2. For instance, the spec¬ 
trum oc |fc| Q with a < 1 in d = 1 requires the inter-site 
hopping oc r~ 1 ~ a which has already been realised in ID 
chains [201121] and 2D arrays[223 of ultracold trapped ions. 
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Utilising fractional a < 1/2 in ID systems also al¬ 
lows one to compare the properties of the unconventional 
disorder-driven transition for \a — 1/21 <C 1 with theo¬ 
retical predictions [3 12] based on the RG approaches con¬ 
trolled by the small parameter e = d — 2a. 

Model. In this paper we focus on a chiral (lacking 
reflection symmetry) ID system described by the Hamil¬ 
tonian 


RG, where l is the (flowing) mean free path. Such be¬ 
haviour of the renormalised disorder strength signifies a 
phase transition between the weak- and strong-disorder 
phases at E = 0. In this paper we verify numerically that 
such a transition persists at all a < 1/2. 

The behaviour of the low-energy density of states near 
a critical point has the generic scaling form (first pro¬ 
posed for 3D Dirac materials in Ref. PI) 


id = |fc|“ sgn k + U(x), (1) 

where k and x are momentum and coordinate and U{x) 
is a short-range-correlated random potential. 

We emphasise that “chiral system” hereinafter means 
a system without reflection symmetry of the quasiparticle 
dispersion ( k —> — k) and should not be confused with the 
concept of a system in a chiral symmetry class[23 ] 124 ] . 

In some sense, such a system is a ID analogue of a 3D 
Weyl semimetal. Indeed, the quasiparticle spectrum con¬ 
sists of two bands: the conduction band (k > 0) and the 
valence band (k < 0), touching at the node k = 0. Quasi¬ 
particles in such a system cannot be localised due to the 
absence of backscattering; the velocity v(k ) = a|fc|“ -1 
never changes sign. In principle, the spectrum of a re¬ 
alistic ID system on a lattice contains an equal number 
of branches with left- and right- movers, due to the con¬ 
tinuity and periodicity of the velocity v(k) (analogously, 
Weyl semimetal has an even number of Weyl points [12]). 
However, for sufficiently smooth disorder elastic scatter¬ 
ing of long-wavelength quasiparticles to states far from 
the node (k = 0) can be neglected and the quasiparticle 
dynamics near the node can be described by the model 

The quenched disorder potential U ( x ) with zero mean 
and a symmetric distribution function preservers the 
E —» — E symmetry of the quasiparticle spectrum and 
the density of states p(E) (with the energy E measured 
from the node), which makes the chiral system partic¬ 
ularly convenient for numerical studies of the disorder- 
driven transition near the node. In contrast, in a system 
with a single band (corresponding, e.g., to the spectrum 
|fc| Q ) quenched disorder generically leads to the renormal¬ 
isation of the band edge and to the formation of Lifshitz 
tails below the band [2], making it hard to define and to 
identify numerically the renormalised edge of the band. 

Disorder-driven transition. For a slightly smaller than 
1/2 the effects of disorder can be analysed using a 
renormalisation-group (RG) approach, controlled by the 
small parameter e = 2a — 1. This RG, previously ap¬ 
plied to Dirac materials [I] [2] [7., [5} I25H25] and to high¬ 
dimensional semiconductors in the orthogonal symmetry 
class [2, repeatedly removes the highest momenta from 
the system, renormalising its properties at lower mo¬ 
menta. Depending on whether or not the character¬ 
istic amplitude W of the random potential exceeds a 
critical value W c , the dimensionless strength of disorder 
7 ~ 1/( k £) flows to larger or smaller values under the 


p(E,W) = Ei 




{W - W c )/E - 


( 2 ) 


with v and z being the correlation-length and dynamical 
critical exponents. In the weak-disorder phase {W < W c ) 
the density of states has the same energy dependency as 
free quasiparticles p(E,W) oc (W c — IU) _1/ (“ _1 )ui _1 , 
while for strong disorder (W > W c ) the density of states 
is smeared and thus energy independent: p(E , W ) oc 
(W — W c )^>. For finite energy E these two regimes 
are separated by a critical region near W = W c with 
p{E) cx E V*- 1 . 

The RG analysis similar to that of Refs. |T] and [2] 
for the model 0 in the one-loop approximation with the 
small parameter e = 2a — 1 yields 

^ = (1 ~ 2a) -1 , (3) 

* = 1/2- (4) 

We emphasise that the dynamical exponent z, Eq. 0, is 
independent of a only in the first order in e and only for 
the ID chiral system under consideration. For arbitrary 
a, not necessarily close to 1/2, the values of the critical 
exponents can be found numerically. 



FIG. 1: (Colour online) The disorder-averaged density of 
states p(E) vs. energy E for a = 0.4 and various disorder am¬ 
plitudes W. For subcritical disorder strength, W < W c « 0.5, 
the density of states vanishes at E = 0, whereas for stronger 
disorder, W > W c , the density of states is finite for all ener¬ 
gies. At the critical disorder strength, W = W c « 0.5, the 
density of states is linear in energy, in agreement with the 
analytical predictions based on one-loop RG calculations. 
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Numerical results. In what follows we present the re¬ 
sults of the numerical simulations that demonstrate the 
existence of the above described disorder-driven transi¬ 
tion for a < 1/2 and its absence for a > 1/2. We also 
analyse the critical behaviour near the transition and ob¬ 
tain numerically the values of the critical exponents. 

To simulate the model ([I]), we use its lattice version 


inter-site hopping element (that we find numerically), for 
long distances \x — x'\ 1 given by the odd power-law 

function [35] 


Jx 


. sgn(a; — x’) T(1 + a) 
\x — x'\ 1+a 27t 


sin 


7T , 

2 (1 + «> 


( 6 ) 


/ ( Jxx'Q^x® J x' “f ^ ^ f ':r ^ J. a ,- , (5) 

x,x f x 

of finite size N with periodic boundary conditions, where 
distances (momenta) are measured in (inverse) lattice 
spacings; J x ^ = J x - x ' = e lk< ' x ~ sc ''>\k\ a sgn k is the 


and U x is the random disorder potential, uncorrelated 
on different sites and described by the Gaussian on¬ 
site distribution with standard deviation W, P(U X ) = 
(W-\/27r) -1 exp[— U%/(2W 2 )]. For the on-site-correlated 
disorder under consideration the ultraviolet momentum 
cutoff[2] K 0 ~ 1 is determined by the lattice spacing. 



FIG. 2: (Colour online) Energy (E) vs. disorder amplitude (W) diagram for the low-energy density of states in a ID chiral 
system with the quasiparticle spectrum = a\k\ a sgnfc. The colour shows the exponent 8 of the density of states p(E) oc E e 
[see the colourbar and Eq. m. For a < 1/2 the density of states has a critical point (E = 0, W = W c ) that separates the 
weak-disorder ( 8 a = a -1 — 1) and strong-disorder (9 = 0) phases at E = 0. The transition disappears for a > 1/2. The grey 
vertical line shows the analytical value of the critical disorder amplitude obtained assuming 0 < 1 — 2a 1. The black crosses 
show the isoline 9=1. The black solid lines show the theoretical crossover energies E = (\W — W c \/W c ) zu between the critical 
region {9 = 1/z — 1) near W « W c and the weak- ( 9 a = a -1 — 1) and strong- (8 = 0) disorder phases, with the exponents v 
and 2 given by Eqs. {3) and Q. 


For each value of a and disorder strength we use the exact diagonalisation method to obtain the spectrum of 
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the system for 100 disorder realisations on a lattice with 
N = 4000 sites and find the parameter 


8{E,W) 


9 In p{E, W) 
9 In E 


(7) 


as a function of the random potential amplitude W and 
energy E, with the results summarised in Fig. [2} 

For sufficiently small a we observe a disorder-driven 
quantum phase transition; the zero-energy density of 
states vanishes RC] for disorder amplitudes smaller than 
a critical value, W < W c , and has a finite value otherwise 
[the dependency p(E) for a = 0.4 is shown in Fig. [I]. 

For sufficiently small a the density of states displays 
three regions with qualitatively different behaviours, that 
touch at a critical point E = 0, W = W c , Fig. [2j (i) for 
low energies and disorder strengths the density of states 
has the energy dependency p(E) oc E 1 / a ~ l of free quasi¬ 
particles with the spectrum oc \k\ a sgn k (0 = 1/a — 1, 
red colour in Fig. [2]); (ii) for low energies and sufficiently 
strong disorder the density of states is constant, 6 = 0, 
(blue colour in Fig. [2]); (iii) near the critical disorder 
strength, W = W c , there is an intermediate critical re¬ 
gion with an almost constant intermediate value of 9. 
These results confirm the existence of the disorder-driven 
phase transition for sufficiently small a. 

In order to accurately verify that the criticality disap¬ 
pears for a > 1/2, we utilise the theoretical predictions 
of the one-loop RG analysis for the critical properties of 
the transition near a = 1/2. Indeed, such RG analysis is 
controlled by the small parameter e = 2a — 1 and thus 
becomes exact when this parameter vanishes. 

The one-loop RG analysis predicts the critical expo¬ 
nents ([3| and Q and the density of states p(E) oc 
E x ! z ~ x = E at the critical disorder strength (W = W c ), 
corresponding to 9 = 1, which can be used to accurately 
identify the critical point. In Fig. [2]the points with 9 = 1 
are shown by black crosses that form a line which con¬ 
tains the critical point and is vertical at low energies. 

The respective value of the critical disorder amplitude 
W c matches well the result (shown by the grey vertical 
solid line in Fig. [5]) 

Wc = [tt( 1 - 2a)/2j* (8) 


of the one-loop RG calculation that we obtain under the 
assumption 0 < 1 — 2a 1 with the ultraviolett mo¬ 
mentum cutoff A'o = 7r (for details of the scheme of the 
perturbative RG calculation see Ref. m ), even for a sig¬ 
nificantly below 0.5. 

These results demonstrate the existence of the criti¬ 
cality of the density of states for a < 0.5 and its disap¬ 
pearance for a > 0.5. For all values of a in the interval 
0.2 ... 0.5 the critical properties of the transition are well 
described by the results of the one-loop perturbative RG 
analysis. 


Conclusion and outlook. In summary, we have demon¬ 
strated that a chiral ID system with the quasiparti¬ 
cle spectrum |fc|“sgnfc with a < 1/2 displays the phe¬ 
nomenology of a high-dimensional disorder-driven phase 
transition. Although all the states in the proposed sys¬ 
tem are delocalised, it exhibits a disorder-driven tran¬ 
sition that manifests itself in the density of states and 
is analogous to the localisation transition near the edge 
of the band of a high-dimensional semiconductor. In 
terms of its symmetries and the critical behaviour of 
observables, the system under consideration presents a 
ID analogue of a 3D Weyl semimetal. The numerical 
values of the critical exponents and the critical disorder 
strength are well described by the results of a one-loop 
perturbative RG calculation. Such a system presents 
a convenient platform for studying high-dimensional lo¬ 
calisation physics and can be used to further investi¬ 
gate strong-disorder conduction in semimetals with delo¬ 
calised states, the interplay of disorder with interactions, 
effects of various disorder symmetries on the transition, 
etc. 

We emphasise, that the unconventional disorder-driven 
transition we studied is not specific to chiral ID systems 
and can be observed (with different critical exponents) 
in any ID chain with sufficiently long-range hopping of 
the excitations, e.g., corresponding to the even disper¬ 
sion oc \k\ a with a < 1/2. However, for systems with 
non-odd dispersions Gc the transition may be harder to 
observe numerically and experimentally due to the renor¬ 
malisation of band edges or nodal points by disorder. 

Implementing the specific ID chiral model in experi¬ 
ments still remains a future research direction. Natural 
candidates are chains of trapped ions, since in those sys¬ 
tems a power-law excitation spectrum oc \k\ a with tun¬ 
able 0 < a < 1.5 has already been demonstrated [201GTl . 
However, ways to generate chiral excitations in these 
chains still have to be investigated. In principle, a non- 
chiral power-law spectrum oc \k\ a is also suitable for the 
observations of the high-dimensional localisation physics, 
but is less convenient for numerical simulations and is 
more sensitive to finite-size effects, that we expect to ob¬ 
scure the disorder-driven transition for non-chiral spec¬ 
tra for small numbers of ions (< 20) used in the current 
experiments. Another candidate for the observation of 
high-dimensional localisation physics is 2D arrays of ions 
in Penning-trap experiments [22j . where a tunable power- 
law spectrum has been demonstrated for about 500 ions 
arranged in a triangular lattice. We leave the analysis of 
finite-size effects and possible realisations of chiral exci¬ 
tations in such systems for future studies. 
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